Casimir potential of a compact object enclosed by a spherical cavity 
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We study the electromagnetic Casimir interaction of a compact object contained inside a closed 
cavity of another compact object. We express the interaction energy in terms of the objects' scatter- 
ing matrices and translation matrices that relate the coordinate systems appropriate to each object. 
When the enclosing object is an otherwise empty metallic spherical shell, much larger than the in- 
ternal object, and the two are sufficiently separated, the Casimir force can be expressed in terms of 
the static electric and magnetic multipole polarizabilities of the internal object, which is analogous 
to the Casimir-Polder result. Although it is not a simple power law, the dependence of the force on 
the separation of the object from the containing sphere is a universal function of its displacement 
from the center of the sphere, independent of other details of the object's electromagnetic response. 
Furthermore, we compute the exact Casimir force between two metallic spheres contained one inside 
the other at arbitrary separations. Finally, we combine our results with earlier work on the Casimir 
force between two spheres to obtain data on the leading order correction to the Proximity Force 
Approximation for two metallic spheres both outside and within one another. 

I. INTRODUCTION 

Casimir forces arise due to vacuum fluctuations of electromagnetic fields in the presence of static or slowly moving 
conductors, or more generally, dielectric or magnetic materials [JQ. The fields obey appropriate boundary conditions 
on the conductors or appropriate constitutive conditions on other electromagnetically active objects, which result in 
induced charges and currents. Due to the quantum nature of the field, the induced charges fluctuate, shifting the 
energy of the vacuum by a finite amount. This difference manifests itself as an interaction — the Casimir force — 
between neutral objects that depends on their sizes, shapes, material properties, and relative orientations. The case 
of perfect conductors is particularly simple: the Casimir force depends only on the geometry of the configuration. 
Analogous Casimir forces can arise from fluctuating scalar or fermion fields in the presence of objects on which 
they obey boundary or constitutive conditions. The electromagnetic Casimir force is a quantum effect observable at 
macroscopic scales. It has been shown to be significant in sub-micron scale devices as well as in the description of the 
interactions of atoms and/or molecules with surfaces, prompting substantial theoretical and experimental investigation 
over the last decade or so. 

In this paper, we report computations of the force on a small polarizable object inside an otherwise empty conducting 
spherical shell as a function of its displacement from the shell's center — the interior analog of the Casimir-Polder 
result. We further give the first exact calculation of the force between a metallic sphere inside a metallic spherical 
shell as a function of their radii and displacement. Although we restrict to metallic surfaces immersed in vacuum/air 
in this paper, the theoretical framework underlying our analysis is universal [2], and our methods have been extended 
to analogous interior configurations involving dielectric objects immersed in dielectric media [3]. Finally, we combine 
our results with earlier work on spheres [4, 5 to obtain first order corrections to the Proximity Force Approximation 
(PFA) for two metallic spheres both outside and within one another. There has been much interest and research in 
computing the Casimir force beyond the PFA [6, 7J; we answer this question in the case of perfectly conducting spheres. 
Some of the results described here were presented in an abbreviated form in an earlier Rapid Communication [8] . 

In the past, there have not been many studies of the Casimir force in closed cavities despite the fact that cavity 
configurations are experimentally realizable. Marachevsky [9 computed the energy of a dilute dielectric sphere and a 
dipole at its center, and recently he studied the interaction of parallel plates inside a cylinder [TO]; Brevik et al. [TT] 
studied concentric dielectric spheres, and Dalvit et al. [12] studied the interaction of a cylinder inside a cylinder. 
Recent theoretical advances [3] (see also precursors Refs. [HEH03]) i n the study of the Casimir force have made it 
possible to analyze a wide variety of geometries and our investigation of the interior case is an example of configurations 
made accessible by these methods. In this paper, we shall only qualitatively summarize the path integral formalism 
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to serve as a reminder for the reader. For an extensive introduction, a review of previous work, and further references 
we refer the reader to Ref. [2]. 

The electromagnetic Casimir energy of an arbitrary configuration C of objects can be calculated from the partition 
function Zc(k), 

m = -^r ^log^K, i 

2tt J q Z {k) 

where Zc(k) is obtained from the Minkowski space functional integral Z(T) = J T>Ae^ s ^ [S[T] is the electromagnetic 
action evaluated from t = to t = T.) after decomposing the fields E into their fourier modes and Wick-rotating 
to imaginary time — k being the imaginary frequency. The restriction n > allows for E(icn) and E*(^ck) to be 
considered independently since E is real. We subtract the Casimir energy at a convenient location (described by 
Zq in the denominator of the log) to remove the cut-off dependent terms in the unrenormalized energy, since such 
contributions arise from the objects individually. 

The spatial configuration of interacting objects manifests itself physically as a continuity condition that the fluctu- 
ating field E obeys on their surfaces. It is possible to trade the constraints on E for fluctuating sources [21 |U [T3] . Then 
the functional integral over the fields becomes free of constraints and can be performed up to a multiplicative constant. 
Since the path integral over the free fields is independent of the location of the objects in space, the multiplicative 
constant is cancelled when dividing by the partition function for the reference configuration. This leaves a functional 
integral over the fluctuating sources on each object a, in which the action is expressed as a functional of the sources 
and the classical field E c i they produce. By superposition, we write E c i = ^ Q E Q!C i, which allows the action to be 
written as a sum over the self- and inter-actions of all the objects in the system. 

With the choice of convenient bases, we can expand E Q c i in terms of multipole fields, generated by the multipole 
moments of the sources induced on S Q . We can express these multipole fields in terms of the transition matrix 
T = (S — 1)/2 of the object under consideration (where S is its scattering matrix) and the multipole sources. The 
self-action of the sources on each object can, therefore, be written entirely in terms of the multipole moments and 
its T-matrix. Similarly, we can express the inter-action of two different objects in terms of their multipoles and a 
translation matrix which relates their coordinate systems in the appropriate bases. Finally, the functional integral 
over the multipole moments of the sources can be performed, leaving an expression for the Casimir energy in terms 
of the objects' T-matrices and the translation matrices. 

We are interested in a situation where one object, the internal (subscript i), is enclosed entirely within a cavity of 
another, the external object (subscript e). The polarizable internal object interacts with fields scattered inside the 
cavity in which it is immersed. Therefore the T-matrix of the external object relevant to the interior case differs 
from the T-matrix that describes scattered waves outside the external object. Following Ref. [2], we denote the 
T-matrices of the internal and external objects by Tf e and T" respectively, where the subscript denotes the object 
and the superscripts denote the relevant scattering amplitude (T ee being the standard T-matrix of an object). For 
conducting boundary conditions on the cavity of the external object, T™ = [TT e ]~ 1 , where T ee would be the standard 
T-matrix for scattering exterior to a conductor in the shape of the cavity. Additionally, the translation matrices in 
the interior problem are different (V instead of U as employed in Refs. [HQS]) because we are interested in quantum 
fluctuations internal to the cavity and external to the internal object. The V-matrices appear because they relate 
regular waves to outgoing waves as opposed to the W-matrices that relate outgoing waves to outgoing waves. With 
these modifications, Eq. ([T]) evaluates to 

£ =* f^ln det ^-^™. (2) 
2ir J dct(Z-T e M T ee ) 

The dielectric properties of the two objects, and the medium separating them inside the cavity are encoded in the 
respective T and V-matrices. The denominator subtracts the energy when the centers of the two objects coincide (as 
opposed to infinitely separated in an exterior problem [4lll6|). In this way we eliminate the cutoff dependent Casimir 
energy of the dipole at the center of the sphere [9]. For a detailed discussion, we refer the reader to Ref. [2]. 

Eq. ||) can be evaluated exactly for certain geometries for which the T and V-matrices are easily calculable in a 
convenient basis. The case of spherically symmetric dielectric objects is particularly simple because their T-matrices 
are diagonal in the basis of spherical wave functions. In this paper, we provide results for a metallic sphere inside 
an otherwise empty metallic spherical shell, while dielectric objects immersed in a dielectric medium are treated in 
Ref. 0. 

The matrix identity IndetM = TrlnM, allows for a simple physical interpretation of Eq. We can express 
the Casimir energy as a series, £ = hc/2-K J dnTr U\f + \M 2 + ...), over the matrix Af = T"Ve,iTi e Vi,e where Af 
describes a wave that travels from one object to the other and back jlj. In general, all terms in this series are 
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important, illustrating the fundamentally non-two-body nature of the Casimir force. The rate of convergence of this 
series depends on the size of the internal object relative to the separation of its surface from that of the cavity. 

First we consider an object that is small compared to the size of the cavity. Then the first term in the series 
expansion of Eq. (|2j), £ — hc/2n J dnTiAf, already gives an excellent approximation to the energy. Furthermore, in 
this limit the Casimir energy is dominated by the lowest frequency contributions from the lowest partial waves in Tf e ■ 
In a spherical basis, the leading terms in the electromagnetic T-matrix are, T^ Vm , ~ k 1+1 +1 and T^ Vm , ~ k 1+1 +2 for 
A 7^ a - , where Z = 1,2, and A and a label the polarizations E (electric) or M (magnetic). The leading contribution 
to the Casimir force comes from the orientation dependent dipole response of the internal object to a dipole field, 

where the internal object can be characterized by its polarizability tensor, ce^f = ^K~ 3 T^/im' ( see S3)- 

Now we fix the external object to be a conducting spherical shell of radius R and define a to be the displacement 
of the center of the internal object from the center of the shell. To leading order in r/R (where r is the typical size 
of the internal object), the Casimir energy can be expressed as [5J, 

^£{a/R) = [f E (a/R) - f E (0)} Tra B + g E (a/R)(2af z - a E x - a E y ) + (E O M) + . . . (3) 

where the polarizability tensor a mm i ~ r 3 has been expressed in a Cartesian basis, and "+..." denotes terms that are 
higher order in r/R. The functions / and g are plotted in Fig. [5] and their functional forms are given in Section [V] 
Eq. ([3]) describes the interaction of a polarizable object (an atom for instance) inside a conducting spherical shell, 
which is analogous to the well-known Casimir-Polder potential [18] . However it differs from the Casimir-Polder result 
in three ways: f,g are non-trivial functions of a/R; the internal object experiences a torque; and the Casimir force 
between the two objects depends on the internal object's orientation. Note that the expansion in Eq. ([3]) is asymptotic 
in r/R at fixed a/R. For a spherically symmetric internal object, the orientation dependent terms in the Casimir 
energy vanish, and corrections to Eq. ^ come from the static quadrupole electric and magnetic polarizabilities, ct^' E 
of the internal object. These corrections are given in the Appendix. 

The opposite extreme from the Casimir-Polder limit occurs when the interior object nearly touches the cavity wall. 
The leading behavior of the Casimir force in this limit is given by the Proximity Force Approximation [19], The PEA 
prediction for the Casimir force between two conducting spheres, whether they are separated or contained one inside 
the other is given by 

^ md ^(d,s,R) = -^^- , (4) 
d^o 360 R + r 

where r and \R\ are the radii of the internal and external spheres respectively and d is the minimum distance between 
their surfaces. By convention we keep r fixed and let R vary. R > corresponds to the exterior problem (the 
spheres are separated); R < corresponds to the interior problem; in the limit R — > 00 we have a sphere opposite a 
plane. The constraint r < \R\ avoids double counting of sphere-sphere configurations. Eq. Q is derived for R > 
by semi-classical methods in Ref. [5D] and the extension to R < is straightforward but its corrections have up to 
now not been known. The planar and exterior problems have been studied in Ref. @] and Ref. [S], respectively (see 
also Refs. [7]). Since most experiments up to now have considered spherical conductors separated by distances 
much smaller than their radii, the first correction in d/r to the PEA is the geometric correction of greatest immediate 
interest. As discussed in Section \FV\ we parameterize the first correction to the PFA by 



Our evaluation of Eq. ^ for conducting spheres has allowed us to combine our results with those of Refs. [4j [5] to 
predict an estimate of the PFA correction coefficient 6\{r/R) appearing in Eq. (p)]) for —1 < r/R < 1. We refer the 
reader to Section llVl for further discussion. 

The rest of this paper is organized as follows: Section [n] provides representations of the vector transition and 
translation matrices relevant to the interior case in a spherical wave basis, followed by an exact computation of the 
Casimir force between metallic spheres in Section [TTTJ In Section [TV] we discuss first order corrections to the Proximity 
Force Approximation for two metallic spheres of arbitrary size based on numerical results in this paper and in [1J |S] . In 
Section |Vj we derive the interior Casimir-Polder result, and study its comparison with the exact results of Section |III| 
Ref. 21 repeats the analysis of this paper for a complex scalar field, which follows by analogy with the vector case. 



II. THEORETICAL BACKGROUND 



As pointed out in the introduction, the Casimir energy of a configuration of compact objects can be evaluated 
in terms of their transition matrices T = (S — T) /2 (where S is the scattering matrix) and translation matrices 
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FIG. 1: Interior Geometry: An object inside the cavity of an external object. We assume that it is possible to choose a 
bounding sphere that does not overlap with the surface of the cavity. The distance between the two origins is denoted by Xi e . 



that relate various coordinate systems relevant to each object. For the interior geometry, we describe the internal 
object by Tf e , the scattering amplitudes for scattering outside its external surface, and the external object by T™, 
the scattering amplitudes for scattering inside the cavity of the external object which contains the internal object. 
We use the subscript M for the medium inside the cavity, and index the dielectric constants e x {icn) and /i x (icK,) by 
x — (i,e,M). Fig. ([TJ is an illustration of the interior geometry. 

Once the geometry and the dielectric properties of the interacting objects are specified, the next step is to calculate 
the T and V-matrices in a convenient basis. We specialize to spherical coordinates and solve the interior problem 
exactly for a sphere contained inside a spherical cavity. Choosing r and R to denote the radii of the internal object 
and the cavity respectively, the T-matrices are diagonal, and are represented in a spherical basis as, 

(j,Mii(nMKr)d r (rii(riiKr)) - fJ.id r (rii(n M Kr))ii(niKr) 

'i.lmMl'm'M — d U'°i 



-> -6u>6. 



HMh{nMnr)d r {rii{ninr)) - Hid r (rki(n M Kr))ii(niKr) 
ii(n M nr) 



ki(n M Kr) 

_ e M ii(n M Kr)d r (rii(niHr)) - e l d r (rii(n M Kr))ii(n i Kr) 

'i.lmEl'm'E ~ "IV "mm' i / \o / \\ a /„; i \\ ■ i \ 

£MKi{nMKr)d r {rii(niKr)) - €id r {rki(nuKr))ii[niKr) 

>-0wd mm r—— (6) 

ci-i-oo o r {rki(nMKr)) 

where n x — ,ye x /i x is the refractive index, and M, E denote the magnetic and electric polarizations respectively. T e n 
can be obtained by the substitutions ii -h- ki, r — > R, and the subscript i — > e for the dielectric constants everywhere 
in Eqs. ([6]). fc; and i; are modified spherical bessel functions of integer index I. 

The spherical wave translation matrices that relate the coordinate systems centered on the internal sphere and the 
external cavity are given by [2 , 



v ie , lWM , lmM = (-ir£w + 1) + ni' + 1) - + 1 )]y T(2f tffw ) +i) +1) 
f 1 1' \" \ f i i' i" \ ,, 

x ( o ) ( - ) *i"(»MK|Xfe|)(-l) Fi„ m _ m ,(X ie ) 



Vie.l'm' E.lmM z 



X- {\f m Bi' m ii m+ x(Ki e ) + A im Bj' m '; m _i(Xie)) , 



V { {HmBl'mnm+xQ 

Vie,l'm'M,lmE = — Vie,l'm'E,lmM ; Vie,l'm' E,lmE = Vie,l'm' M,lmM : (7) 



y— (^ rn Bli m >lm+l{X-i e ) — A im _B;' m '; m _i(X ie )) + zmSp m /; m (Xj e ) 
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where 

iWim(X) = (-l) m J2 VMM + W + 1)(2Z" + 1) (J o o)(m -m' m'- m ) «"(»w«|X|)(-l)''V,» ro _ m ,(*) 

and A^ n = -^/(/ T m)(Z ± m + 1). Xj e is the displacement vector that extends from the center of the internal object 
to the center of the cavity with a = |Xj e |. The translation matrix V e ; is related to Vi e by, 

"-0-0*0-0 (8) 

in the two dimensional space of electric {E) and magnetic (M) polarizations. To simplify calculations, we have aligned 
the z-axis of the cavity along Xj e . The Casimir energy of dielectric spheres immersed in a dielectric medium is derived 
using the above equations in Ref. [5] . We present results for conducting boundary conditions in the following section. 



III. COMPUTATION FOR METALLIC SPHERES 



In this section, we analyze the Casimir interactions of a metallic sphere contained within an otherwise empty metallic 
spherical shell. The Casimir energy is obtained by numerical integration of Eq. ^ using the matrix representations 
given in Section [IT] It depends on the ratio of the radii, r/R, and varies with the displacement a of the centers, 
parameterized by x = a/(R — r). x ranges from zero (the spheres are concentric) to unity (the two spheres touch). 
The Casimir energy varies over many orders of magnitude and indeed diverges as the spheres touch at x = 1. To 
make our graphs easier to read and the subsequent numerical fits easier to perform, we seek to divide the Casimir 
energy by a simple function that captures the growth near the limit x — > 1. This can be accomplished by utilizing the 
Proximity Force Approximation |19j . which accurately predicts that in the limit x 1, the Casimir energy diverges 
as (1 — x)~ 2 . Unfortunately, the leading term in the PFA does not behave correctly as x — > where the force should 
vanish. (The Casimir energy is quadratic in x for small x.) To accommodate both limits, x — > and x — > 1, we 
employ an extension of the PFA which although has no theoretical foundation (beyond the leading term) yet provides 
a definite and convenient function that captures the dominant variation in the exact Casimir energy over the whole 



range of x. We refer to this function as the "full PFA" , and denote it by £{pfa(x)', it is described in detail in Section IV 



(see Eq. (12|). To illustrate our results, we choose r/R = 0.5 and plot R(x) — £ {x) / Sipf\(x) in Fig. 

Given data like that shown in Fig. [2j we obtain the Casimir force by numerically differentiating R(x) with the help 



of the relation, 



-^fPFA(^) -^fPFA^) 

Note that R(x) is determined numerically only up to x = 0.925. Therefore in the range x < 0.9, R'(x) is determined 
by numerical differentiation. For x > 0.9, R'(x) is determined either by differentiating a suitable function that 
extrapolates R{x) or by extrapolating R'{x) itself. Both procedures give identical results (within numerical error, as 
discussed below). The numerical integration and differentiation were performed with MATLAB while all fitting and 
extrapolation procedures were performed with GNUPLOT. 

Although all partial waves contribute to the Casimir energy, partial waves with I < Z max contribute the most; 
where Z max depends on the spheres' relative sizes and separation and grows rapidly as the separation gets small. For 
r/R, = 0.5, the gray curve in Fig. [2] shows the results for ? max < 25. For x > 0.7 the limitation to ? max < 25 is 
inadequate. To obtain results in this range of x it is necessary to include progressively larger values of I. Eventually 
the numerical evaluation of Eq. ^ is limited by our ability to manipulate large matrices. For example Z max > 65 
at x ~ 0.9. To obtain accurate results for x > 0.7 we first compute £ (Z max ) for a sequence of values of Z max . Then 
we fit £ (i max ) to a decaying exponential (which seems to capture the leading behavior at large l max ) of the form 
£{lmax) = £(oo) — ae~^' max , where a and /3 are constants. The resulting function, £(po) is plotted in Fig. [5j It 
smoothly extrapolates to the PFA point at x = 1 as it should. The convergence of £ (l ma x) to £ (oo) is illustrated in 
Fig. [2] (inset). 

At even closer separations, x > 0.925, important contributions to the Casimir energy come from values of I even 
larger than Z f=a 65. It is difficult to evaluate numerically stable values for matrix elements involving modified spherical 
bessel functions ki(x) in the limit x — > with I ~ 65 and above. Therefore our numerical methods are inadequate 
when I grows above 65. But we know that the Casimir energy approaches the PFA limit as x — > 1. Therefore we 
extrapolate the exact data calculated at x < 0.925 to estimate the Casimir energy for 0.925 < x < 1. For example, 
for the case shown in Fig. [2j this is achieved by extrapolating the five data points between 0.825 < x < 0.925 to a 
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FIG. 2: Casimir energy between two conducting spheres. The dark black line shows the Casimir energy, R(x) = £/£{pfa, as a 
function of x = a/(R — r) where a is the displacement of centers. The radius of the inner sphere is fixed at r = 0.57?, where R 
is the radius of the outer sphere. In the limit x — > 1, the Casimir energy approaches the PFA energy, which is marked by the 
dashed horizontal line. £fPFA denotes the 'full' form of the PFA energy discussed in Section [rV] At intermediate separations, the 
Casimir energy is dominated by lower partial waves. For example, the gray line shows that the energy obtained by integrating 
Eq. Q to partial wave order / = 25 is accurate up to x ~ 0.7. The dark line is obtained by extrapolating to I — oo. Inset: 
Convergence at close separations, 0.75 < x < 1. 



function, f(d/r) = l + 9 1 d/r + 9 2 \og(d/r)d 2 /r 2 , where d = R-r-a, X = 1.770±0.034 and 9 2 = 2.272±0.271. Notice 
that 6\ = 6\ — 6*i fPFA where 9\ is defined for all r/R in Eq. ^ and #i,fpFA is calculated from Eq. ( 12 ), and likewise for 
#2 = #2 — #2,fPFA- Therefore, the values of the PFA correction coefficients for r/R = 0.5 can be easily determined. The 
above analysis can be easily repeated for other values of r/R and the coefficients 0\ and 9 2 determined for a range of 
those values. Thus our numerical methods yield sub-leading corrections to the PFA. We perform these computations 
in Section |TV} For more details, the reader is referred to that Section. 

The Casimir force between two conducting spheres depicted in Fig. [3] is calculated by numerical differentiation 
of the exact data points spaced at Ax = 0.025 along the black curve R(x) in Fig. [2] We remind the reader that 
Fig. [2] plots R(x) — £/£fppA as a function of x. The curve shown in Fig. [3] is calculated from R{x) using Eq. (JoJ) . 
The differentiation of R(x) is performed using centered differences for 37 data points between 0.05 and 0.925. For 
x > 0.9, J-/ TfpFA can be determined either by an independent extrapolation or by an algebraic manipulation of the fit 
describing R(x) in that range. We fit a function of the form, h(d/r) — l+6 1 d/2r— 9 2 d 2 /2r 2 — #i,fPFA(^i+#i,fPFA)d 2 /4r 2 
to the four data points between x = 0.825, and x = 0.9 and compare the new extrapolation constants with the ones 
determined in Fig. [2] The function h is determined algebraically from .F/J-ppA, where T is given by Eq. ^ We find 
Q\ = 1.770 ±0.032 and 9 2 — 2.058 ±0.529, which agree with their previously calculated values within their error. This 
achieves two goals: it makes contact with the PFA prediction in Eq. ([5]), and demonstrates that the function f(d/r) 
used to extrapolate R(x) at values of x close to 1 was the correct ansatz for the sub-leading PFA behavior for the 
energy in Fig. [2] 

We have illustrated the numerical evaluation of Eq. ([2| with the case r/R — 0.5. However, the same techniques 
may be applied to determine the Casimir force and energy by numerically integrating Eq. ([2| for all configurations, 
< r/R, < 1. In the following Section, we apply the methods demonstrated in this section to study various r/R con- 
figurations in the limit x — > 1 and determine the PFA correction coefficients 9\ corresponding to those configurations. 
On the other hand, the Casimir force for intermediate values of a/R is studied completely analytically for the range 
of interior configurations r/R — > in Section [v] 
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FIG. 3: Casimir force between two conducting spheres. The black line shows the Casimir force, J- j J-fPFA, between two 
conducting spheres as a function of x — a/(R — r) where a is the displacement of their centers. The radius of the inner sphere 
is fixed at r = 0.5R, where 7? is the radius of the outer sphere. In the limit x — > 1, the Casimir force approaches the PFA. 
J-fPFA denotes the 'full' form of the PFA discussed in Section [TVl 



IV. CORRECTIONS TO THE PFA 



As mentioned in the Introduction, one of the most interesting quantities made accessible by our methods is the first 
non-trivial correction to the Proximity Force Approximation. In this section we extract this correction for the case of 
one sphere within another, and combine it with data from the cases of two separated spheres and a sphere opposite 
a plane, to survey the full range of possible sphere-sphere configurations. 

The leading term in the PFA is given by Eq. Q as discussed in the Introduction. The analytic form of the 
corrections to the PFA is unknown in general, however the case of a sphere facing a plane was treated analytically 
in Ref. j7J. We find that our data can be fitted very well by the first few terms in a power series expansion of the 
Casimir force in d/r, as quoted in the Introduction, 

T - R) " - ssp j£ (' + *' WK) i - ^ & + < 10 > 

as d — > 0. (Remember that r/R < (>)0 corresponds to the interior (exterior) case.) This power series expansion of 
the force requires that the energy include a log(c?/r) term, 

ir 3 frr rU f d d 2 d 2 \ 

m r, R) = - 720^2 ( 1 + e ^/R)- + 2 (r/R) log(d/r)- + a(r/R)- + 0(d 3 /r 3 )j (11) 

where we have adjusted signs in Eq. ( 10 ) so t hat corrections to the PFA energy in Eq. (Ill]) correspond to the 
extrapolation function f(d/r) defined in Section III Note that the term proportional to d 2 /r 7 in Eq. (11) does not 
contribute to the force. 

It is useful to have an estimate, however crude, of the interior Casimir energy over the whole range of d/ R in order 
to scale out the rapid variation that makes it difficult to display £ graphically. To this end we extend the PFA over 
the whole range of d, r, and R. The PFA estimate of £ can be calculated by assuming that each interacting surface 
is assembled out of infinitesimal mirrors spaced at a distance Z(£i,C2) from the other surface, where (Ci : C2) are the 
coordinates of the surface chosen as a convenient reference. This algorithm is ambiguous beyond the leading term in 
1/d because there is no unique way to specify the separation between the surfaces. For definiteness, we extend the PFA 
by taking d to be the distance between the surfaces measured radially outward from the smaller sphere and integrate 
over the surface of the smaller sphere. This can be done for both the interior (r/R < 0)and the exterior (r/R > 0) 
configurations. Note that the restriction y = r/R G [—1,1] covers the full range of sphere-sphere configurations as 




FIG. 4: PFA correction coefficients for spheres. r/R ranges from -1 (interior concentric), to zero (sphere-plane), to +1 (exterior, 
equal radii). The data points correspond to the exact values of 9i, calculated numerically, while the solid black curve is a fit 
(see text). Inset: "interior" and "exterior" geometrical configurations. 



long as r is taken to be the radius of the smaller sphere. The result, which we refer to as the 'full PFA', is given here 
(for compactness) as a definite integral (which can also be obtained in closed form), 

■K 3 hcy 2 f 1 

£lPFA = ~ ~36UR~ / dx I 7 \ 3 + (1 + 7 A3 I . V < 

1 " -' 1 (-(1 + y + yd/r)x + y+^(l + y + yd/r) 2 (x 2 - 1) + l) [i + V) 

ir 3 hcy 2 f 1 dx 
tfPFA = „ nn / r^r, y > 



360i? 



x (v,d/r) f(i + y + y d/r)x -y- + y + yd/r) 2 {x 2 - 1) + 1^ 



where x (y,d/r) = - {l + y \ yd/r)2 - (12) 

(For y < we have subtracted the energy when the spheres are concentric as in Eq. Q.) 

Another option would be to measure the distance d radially inward from the outer sphere. This procedure yields 
a different analytical form for the "full PFA" beyond the leading order, illustrating the ambiguity in defining d. To 
distinguish between the two "full PFA" estimates, we refer to the former as r-based and the latter as i?-based. 



If we expand around d/r = in Eq. (12), we find that 0\ jp-p &.{t / Ft) = —r/R — r/(R + r) — 3, which is continuous 
in the interval r/R e (—1,1]. A similar calculation from the ii-based PFA yields a different, yet continuous, form for 
0i,fPFA, i-e. -{3r/R + r/(r + R) + 1). 

The above discussion suggests that corrections to the PFA extend smoothly from the interior case, — 1 < r/R < 0, 
to the planar case, r/R — 0, to the exterior case, < r/R < 1. Results for the special cases r/R = and r/R = 1 
were presented in Refs. 3] and [5] , respectively. In order to give a full description of the leading correction to the 
PFA, we compute the correction for several additional exterior configurations {r/R > 0) and combine those results 
with the results of Refs. ^ and [S], and with our results for interior configurations in order to obtain a form for the 
PFA corrections over the entire possible range of r/R. Fig. [4] displays our results along with the corresponding "r- 
and -R-based" estimates of 9\ . 

The numerical data in Fig. [4] show a smooth transition from the interior to the exterior configuration. Although 
neither of the "full PFA" estimates describes the data, the r-based PFA has a similar functional form and divergence as 
x — > — 1. Therefore, we fit the data in Fig.|4]to a form motivated by the r-based PFA, 9\{x) = —{kix+k2x/{l+x) + k^) 
and find k\ = 1.05 ± 0.14, fc 2 = 1.08 ± 0.08, and fc 3 = 1.38 ± 0.06. This provides a simple form for the leading PFA 
correction for metallic spheres, one inside the other and both outside, which is relevant for many experiments. Notice 
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FIG. 5: Plot of the functions f M/E (a/R) and g M/E (a/R), denned in Eqs. (14 1 and (151, respectively. 



however, that the actual function 9\(x) is not known analytically and that our fit represents a reasonable choice which 
may not be unique. Our results show that the correction to the PFA has a significant dependence on the ratio of 
curvatures of the two surfaces. The correction is a factor of two larger for two spheres of equal radii than for the 
sphere-plane setup; it vanishes near r/R = —0.5; and it becomes positive and large as r/R — » —1. These effects 
should be taken into account when experimental accuracy has advanced to the point that corrections to the PFA can 
be measured. 



V. CASIMIR-POLDER LIMIT OF THE INTERIOR PROBLEM 



In this section, we derive the Casimir energy of a small polarizable object contained inside a metallic spheri- 
cal shell. Using the matrix identity IndetM = TrlnM we expand the integrand in Eq. Q in a Taylor series, 
£ = hc/2n J dnTy (j\f + ^J\f 2 + ...), in N — T"V et iTf e Vi, e where each matrix multiplication by N describes the prop- 
agation and reflection of a virtual photon from the pair of objects [4]. In the exterior case, both objects can be taken 
to be small compared to the separation between them, and can be approximated by their static dipole polarizabilitics 
in this limit, leading to the famous Casimir-Polder force between polarizable molecules. In the interior case studied 
here, the enclosing cavity must be larger than the interior object, and cannot be approximated by its lowest frequency 
electromagnetic response. Therefore, any asymptotic expansion of Eq. |2| requires that the shape and material prop- 
erties of the enclosing cavity be specified and its T-matrix be calculable. We obtain a useful expansion by requiring 
the size of the internal object to be small enough that it be adequately described by its electric and magnetic dipole 
polarizability tensors. Thus the most general interior configuration to which an analytic expansion applies is that of a 
small dielectric object immersed in a dielectric medium inside a dielectric cavity for which the scattering amplitudes 
are known analytically in some partial wave basis. Our results on the Casimir-Polder limit of a small internal object 
inside a metallic spherical shell were reported in Ref. [5], and stable three dimensional configurations of objects in 
spheroidal drops of liquid or metallic shells were presented in Ref. [3] . 

Here we consider a metallic spherical cavity of radius R enclosing a small object of typical linear dimension r, and 
use the T-matrix representations in Eqs. ([6]). Substituting the lowest order approximation to the interior object's 
T-matrix, 7~ ee , but keeping all partial waves in the T-matrix of the exterior sphere, T™, we find the first term in the 
expansion in r/R, with coefficients that are non-trivial functions of a/R. This expansion was already quoted in the 
Introduction, 

^£(a/R) = [f E (a/R) f E (0)} Tra £ + g E (a/R)(2af z - a E x - a E y ) + (E O M) + . . . (13) 



where a denotes the displacement of the internal object along the z-axis from the center of the cavity and a denotes 
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0.25 



FIG. 6: Comparison of the interior Casimir-Polder result with the exact Casimir energy predicted by Eq. Q for conducting 
spheres. Plotted along the y-axis is the fractional error (£ — £cp)/£ as a percentage, where £cp is calculated from Eq. (181. 
For each value of the separation a/R (denoted by point markers listed above), we have plotted the fractional difference between 
£ and £cp at the leading order 0(r 3 /R 3 ) (gray) and next-to-leading order 0(r 5 /7? 5 ) (dark black) as a function of the internal 
sphere radius r/R. 



the dipole polarizability tensor expressed in a Cartesian basis. The coefficient functions f p and g p , plotted in Fig. [HJ 
can be expressed in terms of modified spherical Bessel functions i v and fc„ as T 



f E (0 



dxx 3 V" 
1=1 



( E (x) 



2<r2 



(a + i)^_i(^) + «m(^))-CrW 



2(21 + 1 



(14) 



(0 = jT dxx 3 2^T) (^V^i(^) + - 31(1 + l) lH (xO I1+1 (xO) 



(15) 



and / M and g M are obtained by substituting (i? <-> M) in the above equations. The functions (^l E are given by 



h(x) 



h(x) + xk[(x) 
ii(x) + xi'i(x) 



(16) 



An analogous result — a series in r/R with coefficients that are non-trivial functions of a/R — can be obtained 
for an analogous "exterior" configuration: a polarizable object facing a metallic sphere of radius R, see Ref. [8 . Both 
results differ from the classic Casimir-Polder result, which describes an object facing a conducting plane, in three ways: 
f,g are non-trivial functions of a/R; the polarizable object experiences a torque; and the Casimir force between the 
two objects at leading order depends on the polarizable object's orientation. Specific features of the above orientation 
dependence are explored by specifying the internal object to be a dielectric spheroid in Ref. j8]. For example, it is 
demonstrated that the orientation dependence, like the leading PFA correction, has a smooth continuation from the 
interior to the exterior configurations; a 'cigar-shaped' spheroid prefers to align itself perpendicular to its displacement 
vector from the center of the shell; and a 'pancake-shaped' spheroid prefers to align its two large axes perpendicular 
to its displacement vector. 

For a spherically symmetric internal object, the first non-trivial correction, which is 0(r 5 /R 5 ), to Eq. (13) can be 
easily evaluated. This is possible because the T-matrix for the spherically symmetric object is diagonal in a spherical 



11 



basis and independent of the azimuthal index, m, T,, 



Act 

ml'm' 



5w8 mm '8 Xcr T{' ■ For each polarization and I, the leading 



term in the T- matrix is proportional to a multipole polarizability, af ~ r 2l+1 , which characterizes the low frequency 
electromagnetic response of the internal object, and which can be computed for simple geometries or measured for 
any compact conductor 0], 



77 



21 



(-ly-^ + iK , 

l(2l + l)U(2l- 1)!! 



0(k 3 



for a = E or M. Substituting Eqs. ^ and Q in Tr N we find the Casimir energy up to 0(r 5 /R 5 ) to be, 



2nR r h "( /m a i 



+ /if (a/R)-4r + (Af o £) + 0(r7i? b ) 
it 



(17) 



(18) 



where h^^ E are proportional to f M / E (a/ R) defined in Eq. (14) up to a numerical factor. The exact functional forms 



of the coefficient functions hi and hi are given in the Appendix. Each can be expanded as a power series in a 2 / R 2 
which converges for a 2 /R 2 < 1. The Casimir force can be calculated by differentiating hi, hi with respect to a. 

To examine the usefulness of the expansion in Eq. (18), we compare its predictions with the exact numerical result 
for £ following from Eq. pi) for the case of an internal metallic sphere of radius r, for which af" 1 — —lr 2l+1 /(l + 1) 



and 



„22 + l 



Fig. 



lst-order data (in blac 



6] plots the fractional errors A£ = (£ — £cp)/£ vs. r/R for various separations a/R. The 
;k markers) include contributions from Q( R~ 3 ) terms only, while the 2nd-order data (in gray 
markers) include coefficients from h2(a/R) at 0(R~ 5 ) in Eq. (18). Many trends are visible in this graph: For example 
the interior Casimir-Polder result through second order is accurate to more than 99% of the exact answer for all 
r/R < 0.1 for < a/R < 0.4. Another feature worth noting is that for a given value of r/R, \im a / R ^ Q A£ is not 
zero. Both the exact Casimir energy, £, and the interior Casimir-Polder approximation, £qp, vanish like o 2 as a — > 0. 
(Remember the value of each at a = has been subtracted.) Notice however, that the Casimir Polder approximation 
is an expansion in r/R, not a/R, and therefore each term in the expansion, Eq. (18), contributes, albeit with smaller 
magnitude, at a = 0. Fig. [6] shows that the interior Casimir-Polder expansion becomes exact for an arbitrarily small 
polarizable object (an atom or a molecule) inside a conducting spherical shell. 



VI. CONCLUSIONS 



We have studied the electromagnetic Casimir problem for a compact object contained inside a closed cavity of 
another compact object. Using the scattering formalism, we express the Casimir energy between the two objects in 
terms of their T-matrices, and translation matrices that relate the coordinate systems appropriate to each object. 
Then we specialize to the case when both objects are conducting spheres, and illustrate our methods and results by 
evaluating the Casimir energy for the case r/R = 0.5. The Casimir force for this sphere configuration is calculated 
by numerically differentiating the energy with respect to the spheres' separation. 

We have also calculated the analog of the Casimir-Polder expansion for an object contained inside a metallic spherical 
shell. The Casimir energy can be expanded as an asymptotic series in r/R (the leading term being proportional to 
r 3 /i? 3 ), where R is the radius of the spherical cavity, and r is a length characterizing the size of the internal object. 
There are certain novel features of this result: the coefficients are non-trivial functions of a/R (where a denotes the 
separation of the center of the internal object from the center of the shell) that is, they are represented by infinite sums 
of modified spherical bessel functions; the Casimir force at leading order depends on the orientation of the internal 
object; and the internal object experiences a torque. Additionally, we have calculated the coefficient functions for the 
leading two terms for the case when the inner object has spherical symmetry. A comparison of the 'interior Casimir- 
Polder' expansion (up to 0(R~ 5 )) with the exact energy (calculated numerically) for various sphere configurations 
shows it to be accurate to more than 99% of the exact answer for a/R < 0.4 through r/R < 0.1 where a is the 
displacement of their centers. 

The methods demonstrated in this paper can be applied very easily to calculate the Casimir force between dielectric 
spheres at all separations, although they become computationally intensive when the spheres are nearly touching. In 
this limit, we make contact with the Proximity Force Approximation (PFA). A careful examination of the approach 
to the PFA allows us to calculate leading order corrections to the PFA limit. We then combine our studies of various 
interior sphere configurations at close separations with previous work on metallic spheres exterior to each other and 
a conducting sphere facing a mirror to analyze the leading order PFA corrections for all sphere-sphere configurations. 

This work was supported by the NSF through grant DMR-08-03315 (SJR), DFG through grant EM70/3 (TE), 
the U. S. Department of Energy (D.O.E.) under cooperative research agreement #DF-FC02-94ER40818 (RLJ), and 
MIT's undergraduate research opportunities program (UROP) (SZ). We thank Noah Graham and Mehran Kardar 
for many useful conversations and suggestions regarding this work. 
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Appendix A: Interior Casimir-Polder coefficient functions 



We give the coefficient functions hf 1 
modified spherical bessel function i\. hf and 



and hf 



that appear in Eq. ( 18 ) in terms of and £f defined in Eq. ( 16 ) and 



h% can be obtained by substituting M E in the following equations. 



h?(0 



dxx 




2 e 



21 + 1 



(n-xixO-ii+iixOYt-^ix) 



(Al) 



^(0 



dxx 5 ( ^ 

x*e 



i=i 



»m, n (* " + l)(2i + 3)i?_ 2 (atf) + i(* + 2)(21 - l)zf +2 K) + (3/ + 3/2)^(^) 



3(2? + 1 
+ (i-l)(i + 2) 



1 / 1 — Z 



4 V 2/ - 1 



ii-2(a£) 



2Z- 1 



6(4/(i + l)-3) 
Z + 2 



2(2/ — 1) 



41{l + l) -3 
2/ + 1 



k{x£) 



41(1 + 1) -3 



»i(a:f) + 



2? + 3 
1 

2(2^ + 3) 



iz+2(^) 



-gCf (x) 



(A2) 
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